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I analyse the identified single-particle pt spectra and two-pion Bose-Einstein correlations 
^ I from RHIC. They indicate a massive transverse expansion and rather short lifetime of the 
system. The quantitative analysis in framework of the blast-wave model yields, however, 
^ ■ unphysical results and suggests that the model may not be applicable in description of 
•jH . two-particle correlations. I then discuss generalisations of the blast-wave model to non- 
central collisions and the question how spatial asymmetry can be disentangled from fiow 

0^ ' asymmetry in measurements of V2 and azimuthally sensitive HBT radii. 
O ' 

P : 1. INTRODUCTION 

When studying nuclear collisions at highest energies we are interested in properties of 
' strongly interacting matter. A clear signal of existence of an extended piece of matter 
^ ■ is its collective expansion which is not present in simple nucleon-nucleon collisions. This 
can be deduced from the slopes of identified hadronic single-particle pt spectra: while the 
slope is universal for all particle species in proton-proton collisions, the spectra become 
^ I fiatter with increasing particle mass in collisions of heavy ions. 

' Hadrons interact strongly an so they can only decouple from the system when it becomes 
dilute enough. Their spectra are fixed at the moment of freeze-out and carry information 
about the phase-space distribution in the final state of the fireball. We want to reconstruct 
this information and get an idea about the collective evolution of the fireball by comparing 
its final state to its initial state which is known from the energy and the collision geometry. 

The process of freeze-out is very complicated Here I will simplify it and assume 
that all particles decouple along a specified three-dimensional freeze-out hypersurface j2] . 
The phase-space distribution at the freeze-out hypersurface will be parametrised by the 
so-called blast-wave model. 

I will first introduce the model and analyse single-particle pt spectra and two-pion 
correlations from RHIC. Then I show how the model is generalised to non-central collisions 
and focus on calculation of f 2 and azimuthal dependence of HBT radii. 

2. THE BLAST- WAVE MODEL FOR CENTRAL COLLISIONS 

I will assume that in the end of its evolution the fireball is in a state of local thermal 
equilibrium given by temperature T and chemical potentials fii for every species i. The 
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freeze-out hypersurface will stretch along the r = — t? = tq = const hyperbola^ 
and will not depend on the transverse radial coordinate r. There is longitudinally boost- 
invariant expansion, i.e. Vz = z/t. The velocity field 

= (coshp cosh 7], cosip sinhp, smi/j sinh p, cosh p sinh?]) . (1) 

Here, p stands for the transverse rapidity which will be assumed to depend linearly on r 

p = Por/R. (2) 

Note that this linear dependence was found for small transverse velocities also in hydro- 
dynamic simulations jH]. In the last equation, R is the transverse size of the fireball. 
The transverse density distribution is uniform for r < R. The (Wigner) phase-space 
distribution of particles at freeze-out is then characterised by the emission function^ 



X ,4 COsh(?/ - T]) 

S{x,p)dx = -— 



exp ( ) ± 1 



e{R-r) 



where the term rrtt cosh(?/ — r/) comes from the Cooper-Frye pre-factor da^p^ f2j and the 
exponential in r represents smearing in the freeze-out time. The energy in the statistical 
distribution is taken in the rest-frame of the fiuid: E = p ■ u. This introduces "coupling" 
between momentum of the emitted particle and the fiow velocity of the piece of fireball 
where it comes from. The strength of this coupling is controlled by 1/T. 

2.1. Calculation of spectra 

The single-particle spectrum is obtained from the correlation function as 

d^N C . 

E^= / d''xS{x,p). (4) 



dp^ J 

For low pt, the inverse slope^ of the spectrum is roughly given by j3] 

n = T + m{vtf. (5) 

Here, (vf) is the mean transverse collective expansion velocity. Thus the freeze-out tem- 
perature and the transverse expansion cannot be both determined from a measurement 
of only one single-particle spectrum. However, they can be determined from identified 
spectra of different species. 



"'^As space-time coordinates I will be using longitudinal proper time r = Vt^ — (z is the Cartesian 
coordinate in beam direction) and space-time rapidity 77 = l/21n((t + z)/(t — z)). In the plane transverse 
to the beam, radial coordinates r and are used. 

^Momentum is parametrised with the help of rapidity y, transverse momentum pt, transverse mass 
nit = \/ rn? -\- pf, and the azimuthal angle </> such that — (mjcoshy, pt cos cj), pt sin (j), mts'mhy). 
•^The inverse slope appears as a parameter of a fit to the rrit spectrum with the function 
Af exp{—mt/T^,). 
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The full expression for the single-particle spectrum reads [5] 

dro f mtCOsh{p{r))\ /pt sinh(p(r)) 
T J^^^ T 

Note that in the model given by emission function eq. Q, G{r) = 9{R — r) and the second 
term vanishes because dT^/dr = 0. Here I display a more general expression because I 
want to discuss all effects which influence the single-particle spectrum: 

1. Freeze-out temperature and the transverse expansion velocity. 

2. Contributions to particle production from resonance decays. This is not included 
in the formula © but is important in particular for pion production (and will 
be included in calculations below). In this connection also chemical potentials of 
the individual resonance species are important as they determine the number of 
resonances and thus the fraction of pions which stem from their decays. 

3. The shape of freeze-out hypersurface enters through the term dro /dr. 

4. Also "finer details" like the profile G{r), spatial dependence of the temperature, the 
radial profile of the transverse velocity, etc. can influence the spectrum |4|tij . 

It is important to realise that two models which differ in points may not give the 
same T and (vt) when fitted to single-particle spectra. For example, in the single-freeze- 
out model [ZI To = V const + z^ + r^ and the second term on the r.h.s. of eq. (jHI) gives a 
contribution. Thus the temperature obtained from fits with that model cannot be directly 
compared to the one which will obtained below |H]. 

2.2. Calculation of HBT radii 

I will also calculate the HBT^ radii in Bertsch-Pratt parametrisation (see e.g. [H] for 
more details on HBT interferometry) . They appear as width parameters of a Gaussian 
parametrisation of the correlation function 

Ciq, - 1 = A exp {-RliK)ql - R',iK)q', - RKK)qf) , (7) 

where q and K are momentum difference and the average momentum of the pair, respec- 
tively, and the i?'s are the HBT radii. They measure the lengths of homogeneity ^Uj, i-e., 
sizes of the part of the source which produces pions with specified momentum. The sizes 
are measured in three directions: longitudinal is parallel to the beam, outward parallel 
to transverse component of K, and sideward is the remaining Cartesian direction. The 
parametrisation (|7j) is valid in the CMS frame at mid-rapidity of central collisions; other- 
wise terms mixing two components of q may appear (like Rliqoqi, for example; see P for 

^The acronym HBT comes from names of the scientists who apphed the method first in radioastronomy 
in the 1950s: R. Hanbury Brown and R.Q. Twiss. 
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more details). The HBT radii will be determined as 

Rl{K) = ((i-At)'), (8a) 

R'siK) = (f), (8b) 

RiiK) = {z'). (8c) 

In eqs. dH)) x, y, z, t stand for the Cartesian space-time coordinates, z is the longitudinal 
coordinate and x the outward coordinate. Averaging and the tilde are defined as 

Because the lengths of homogeneity depend on momentum, the HBT radii show depen- 
dence on Kt. For R'^{Kt), an approximate formula^ says jl] 
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^'= 1^ M/T ^^^^ M, = ^m^ + K?. (10) 

Thus the temperature T and transverse expansion measured by po are coupled together 
here. However, recalling eq. (0), they can be disentangled by analysing both the HBT 
radii and the single particle rrit spectrum. 

3. THE FREEZE-OUT STATE IN CENTRAL COLLISIONS 

I will fit the model to single-particle spectra and HBT radii measured at RHIC. First, 
each of the identified spectra of pions, kaons, and protons of both charges will be fitted 
individually. This provides a consistency check for the assumption that all these species 
freeze-out simultaneously: if the results of the fits are incompatible, this assumption 
is invalid. After that, I will fit the HBT radii and look again whether they can be 
accommodated with the same model as the single-particle spectra. 

Resonance production of pions will be taken into account for spectra but not for HBT 
radii. A study of the infiuence of resonance production on HBT radii indicated that in 
the presence of transverse fiow they are changed only marginally |llj . 

Chemical potentials of the resonance species are determined in accord with the model 
of partially chemically frozen gas [12] . It assumes that after the chemical freeze-out [13] 
the effective numbers of particles decaying weakly (and thus slowly) stay fixed while the 
strong (and therefore fast) interactions stay in equilibrium. Effective number of any given 
species includes this species plus particles which can be produced on average from decays 
of all present resonances. Adjectives "fast" and "slow" refer to comparison with the 
typical time scale of the fireball evolution. As an example: effective number of pions 
= Nt, + 2Np + A^A + • • • , and chemical potential of A will be /ia = /Utt + /^at because 
it is in equilibrium with pions and nucleons. 

Au-|-Au at ^/s = ISOAGeV. When fitting the various single-particle spectra ^], no 
overlap of the fit results is found at la level. The contours in freeze-out temperature 
T and the average transverse expansion velocity (vt) are overlapping at 95% confidence 
level at T ^ 120 - 140 MeV and {vt) ^ 0.45 - 0.5. 

^Actually, this formula was derived in a slightly different model, where Rg was the transverse size. 
Nevertheless, it can be used for the qualitative arguments here. 
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Figure 1. Two fits to the HBT radii measured in Au+Au collisions at ^/s = ISOAGeV. 
Dashed lines show a fit with values that also reproduce the single-particle spectra: T = 
121 MeV, {vt) = 0.51, R = 12.1 fm, tq = 4.3fm/c, At = 4.4fm/c. Solid lines show the 
best fit to 7r+7r+ correlations: T = 33 MeV, (vt) = 0.73, R = 24.1 fm, tq = 21.3 fm/c, 
At = 1.1 fm/c. 



Can we fit the HBT radii with these values? In Figure ^ I plot the HBT radii measured 
by STAR J3] and PHENIX ^Hl collaborations together with two fits. It turns out that 
Rs{Kt) is so steep that — within this model — it rules out the temperatures obtained from 
spectra analysis. A lower temperature is favoured, the best fit being obtained for T = 
33 MeV! But this is not the freeze-out temperature! The unphysically low value of the 
parameter T just indicates that in order to reproduce the observed strong Kt dependence 
of HBT radii the coupling of momentum to the local flow velocity in the emission function 
eq. © must be very strong. 

Au-|- Au at ^/s = 200 AGeV. Results from fits to identified single-particle spectra of 
pions, kaons and protons |17|18j overlap at la level at T ^ 115 MeV and (vt) ~ 0.57. 
The fit to HBT radii (Fig. ^ with these parameters is marginally good. The best fit 
is again obtained with an unphysically low freeze-out temperature. The distinguishing 
power between the two models will be improved when final data will be fitted. Note that 
a fit to final results from PHENIX seem to be in better agreement with the conventional 
parameter values j21j . 

A very low freeze-out temperature is also obtained when the complete collection of 
HBT radii form Pb+Pb collisions at projectile energy of 158 AGeV is fitted. That fit is 
not shown here due to lack of space. 

Conclusions from central collisions. The parameters obtained in the fits to Kt 
dependence of HBT radii have no direct physics interpretation. The low temperature 
just indicates the need for very strong coupling between momentum and flow within the 
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Figure 2. Fits to the HBT radii measured in Au+Au collisions at = 200AGeV. 
Preliminary data are read off the figures presented at the Quark Matter 2002 conference 
|19|20|j . Dashed lines show a fit with values that also reproduce the single-particle spectra: 
T = 115 MeV, {vt) = 0.57, R = 13.3 fm, tq = 4.7fm/c, At = 4.7fm/c. Solid lines show 
the best fit to the HBT radii: T = 28 MeV, {vt) = 0.71, R = 27.7 fm, tq = 24.9 fm/c. 
At = 1.7 fm/c. 



blast-wave model. The model seems not to describe the realised scenario of freeze-out. 

In spite of all these failures, one has to notice that the measured Rs, which is the 
rms radius of the source, indicates a strong transverse expansion since the original rms 
radius of Au nucleus is about 3.5 fm. If we model the nucleus and the fireball by a box- 
like transverse distribution, this is an expansion from a radius of 7fm to about 12 fm 
(in the "conventional" models with T ^ 100 — 120 MeV). The total lifetime in Bjorken 
longitudinal expansion scenario can be estimated from tq + At ^ 9 fm/c. It remains 
puzzling how the fireball expands transversely so much in such a short time without any 
original transverse flow. 

4. NON-CENTRAL COLLISIONS 

Next I generalise the model for non-central collisions. The transverse profile becomes 
ellipsoidal instead of circular, so this changes the 6'-function from eq. Q into 




e{R-r) \R-Rj— + j-] , R, = aR, Ry = R/a. (11) 



Here I define the direction of the x-coordinate to be parallel to the impact parame- 
ter, and y-coordinate as perpendicular to the reaction plane. The spatial anisotropy is 
parametrised by the parameter a. The transverse flow rapidity — given by eq. © for 
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central collisions — will also vary with the angle ip with respect to the reaction plane 
P(^)V') = ^J cos2 ip + a"^ sin^ ip po (1 + p2 cos(2?/;flow)) • (12) 

I will discuss two models which differ in the azimuthal dependence of p{r,ip) |22j . 
Model 1: Transverse velocity is perpendicular to the surface and V'flow = Arctan| j23j . 
Model 2: Transverse velocity points outwards in radial direction, ipaow = '4'- 

4.1. Elliptic flow coefficient t>2 

Elliptic flow can be generated by two different mechanisms in our model. Firstly, the 
fireball can be spatially ellipsoidal, but have the same transverse velocity at the outer 
shell. Secondly, the fireball can have spatially circular cross-section, but the expansion 
velocity can be higher in one direction than in the other. It would be interesting to 
disentangle these two effects from the data. For example, if we can learn about the shape 
of the fireball, we obtain an indirect information about the lifetime because the fireball 
starts as out-of-plane elongated and via expansion this shape becomes more and more 
circular. An out-of-plane elongated (i.e. a < 1) fireball at freeze-out is thus an indication 
of short lifetime. 

The coefficient f 2 can be calculated in both used models; it turns out that the results 
obtained for Model 1 map exactly to those for Model 2 if one replaces a — > 1/a |22j . 
Hence, without the knowledge about which model applies better for the description of the 
freeze-out, one cannot say anything about whether the fireball is in- plane or out-of-plane 
elongated at freeze-out just by looking at V2. If the model is fixed, V2 depends on a and 
P2 in a different way for different species, so one should be able to disentangle them from 
a measurement of V2{pt) of more than one identified species. 

4.2. Azimuthally sensitive HBT 

Since the HBT radii are measured along directions specified by the orientation of the 
transverse momentum, they depend explicitly on the angle between the transverse mo- 
mentum and the reaction plane. In addition to that, homogeneity regions depend on the 
momentum of produced particles and this introduces an implicit azimuthal dependence. 
(See jD] for more on azimuthal dependence of HBT radii.) 

One quantifies the azimuthal dependence of HBT radii by decomposing them into 
Fourier series. Up to leading oscillating terms this reads 

R'oi^) = Ko + 2^0,2 cos(20) + . . . , R^icjy) = Rl, + 2Rl, cos(20) + . . . , 

i?L(0) = 2i?L,2 sm(20) + . . . , Rf{<p) = Rl, + 2Rl, cos(20) + . . . . ^'"^ 

We will be interested in the dependence of second-order amplitude of the outward radius 
i?o2 the sideward radius i?^2 the spatial anisotropy a and flow anisotropy p2. We 
would like to eliminate their scaling with the absolute transverse size of the fireball, so 
we will divide them by the respective average terms R^q- 

This is plotted in Figure El The good news from this figure is that almost in all cases 
the normalised oscillation amplitude Rj2l^o depends mainly on spatial anisotropy and 
only weakly on flow anisotropy. The only found exception is i?^ in Model 2: its oscillation 
is dictated by flow asymmetry at high pt. At low pt, the homogeneity regions are mainly 
characterised by the shape of the fireball. At higher , they become more sensitive to the 
flow. In case of Model 2 this can even lead to a change of sign of the amplitude i?^ 2 E3 • 
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Figure 3. The second order Fourier terms normalised by the average HBT radii 1^2/^1,0 
(see eq. (HHI)) as a function of a and p2- Upper row: outward radius; lower row: sideward 
radius. Thick lines show where there is no second-order oscillation of the radius as a func- 
tion of (j). Consecutive lines correspond to increments by 0.1. Left panel: calculation with 
Model latpt = 300MeV/c (left) and pt = 900MeV/c (right). Right panel: calculation 
with Model 2 at = 900MeV/c; results at pt = 300MeV/c are similar to Model 1. 



4.3. Comparison with data 

By comparing to data one can try to distinguish whether we observe an in-plane or out- 
of-plane elongated fireball and which of the two introduced models is better. In Figure ID 

1 show V2{pt) compared to curves calculated with Model 1. Within this model we obtain 
an out-of-plane elongated source (a < 1), i.e., a source which remembers its original 
deformation. Recall, however, that the same theoretical curves can be obtained in Model 

2 if we change a into 1/a, so we would have an in-plane elongated source. 

The key to resolve this ambiguity is a fit to the data on azimuthal dependence of HBT 
radii. I can fit their oscillation well with Model 1 if I assume parameters from the fit 
to V2', this is shown in Figured I checked that in Model 2 the amplitude of oscillations 
is opposite to what the data show [22] • The conclusion is that the observed fireball is 
out-of-plane elongated and that Model 1 reproduces data better than Model 2. 

With the used "conventional" values of temperature and average radial fiow gradient po 
I could not reproduce the absolute size and Kt dependence of the sideward HBT radius, 
but note that this is the same kind of problems as was observed in central collisions. 
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Figure 4. Elliptic flow coefficient 
f 2 measured by the PHENIX col- 
laboration in Au+Au collisions 
at = 200 AGeV for the cen- 
trality class 20-30% |21 as a 
function of pt for a pions and 
kaons (squares) and for identi- 
fied protons (circles). Theoreti- 
cal curves are obtained with the 



0.5 1 1.5 2 Model 1 using parameter values 

Pt [GeV] indicated in the figure. 



5. CONCLUSIONS 

We do see a strong collective transverse expansion at RHIC (and at SPS) because the 
transverse size of the source is much larger than the original nuclei. The total lifetime 
seems to be short: this is indicated by the size of Ri when interpreted in Bjorken scenario 
and by the fact that in non-central collisions we observe an out-of-plane elongated source 
so there wasn't enough time to wash away this anisotropy. On the other hand, the lifetime 
must also be long enough to allow the fireball for a massive transverse expansion. This, 
together with the failure of the blast-wave model in fitting the HBT radii, indicates that 
the model does not provide a satisfactory description of the decoupling of particles. As 
the freeze-out temperature is always obtained in a specific model there is no conclusion 
about its particular value. 



Acknowledgements 

I thank the organisers for the invitation and for the stimulating environment created 
at the conference. I am grateful to Scott Pratt and Evgeni Kolomeitsev for inspiring 
discussions and to Mike Lisa and Fabrice Retiere for critical reading and comments to the 
manuscript. This research was supported by a Marie Curie Intra- European Fellowship 
within the 6th European Community Framework Programme. 



REFERENCES 

1. V. Magas, these proceedings. 

2. F. Cooper and G. Frye, P hys. Rev. D 10 (1974) 186. 

3. P. F. Kolb, |arXiv:nucl- th/0304036| 

4. T. Csorgo and B. Lorstad, Phys. Rev. C 54 (1996) 1390. 

5. E. Schnedermann, J. Sollfrank and U. Heinz, Phys. Rev. C 48, 2462 (1993). 

6. B. Tomasik and U. Heinz, Eur. Phys. J. C 4 (1998) 327. 

7. W. Broniowski and W. Florkowski, Phys. Rev. Lett. 87 (2001) 272302. 



10 




Figure 5. Azimuthal dependence of the HBT radii at midrapidity from Au+Au collisions 
at 200 AGeV and centrality class 20-30% measured by the STAR Collaboration |25j . 
From highest till lowest the data are taken for Kt values: 0.2, 0.3, 0.4, and 0.52 GeV/c. 
Theoretical curves are calculated with Model 1 with the parameter values: T = 120 MeV, 
po = 0.99, p2 = 0.035, a = 0.94646, R = 9.4 fm, tq = 5.0fm/c, At = 2.9fm/c. 



8. U. Heinz, talk at the RIKEN BNL Workshop, Nov. 17-19, 2003, 
|Etp://t onic.physics.sunysb.edu/flow03/ 

9. B. Tomasik and U.A. Wiedemann, in Quark-Gluon Plasma 3, R.C Hwa and X.- 
N. Wang eds.. World Scientific, 2003, gepj)h/0210250 

10. A.N. Makhlin and Y.M. Sinyukov, Z. Phys. C 39 (1988) 69. 

11. U.A. Wiedemann and U.W. Heinz, Phys. Rev. C 56 (1997) 3265. 

12. H. Bebie, P. Gerber, J.L. Goity and H. Leutwyler, Nucl. Phys. B 378 (1992) 95. 

13. P. Braun-Munzinger, D. Magestro, K. Redlich and J. Stachel, Phys. Lett. B 518 
(2001) 41. 

14. K. Adcox et al. [PHENIX Collaboration], Phys. Rev. Lett. 88 (2002) 242301. 

15. C. Adler et al. [STAR Collaboration], Phys. Rev. Lett. 87 (2001) 082301. 

16. K. Adcox et al. [PHENIX Collaboration], Phys. Rev. Lett. 88 (2002) 192302. 

17. J. Adams et al [STAR Collaboration], Phys. Rev. Lett. 92 (2004) 11230. 

18. S.S. Adler et al. [PHENIX Collaboration], Phys. Rev. C 69 (2004) 034909. 

19. M. Lopez Noriega for the STAR Collaboration, Nucl. Phys. A715 (2003) 623c. 

20. A. Enokizono for the PHENIX Collaboration, Nucl. Phys. A715 (2003) 595c. 

21. F. Retiere, J Phys. G 30 (2004) S827. 

22. B. Tomasik, |arXiv:nucl-th/0409074, 

23. F. Retiere and M.A. Lisa, |arXiv:nucl-th /0312024, 

24. S.S. Adler et al. [PHENIX Collaboration], Phys. Rev. Lett. 91 (2003) 182301. 

25. J. Adams et al [STAR Collaboration], Phys. Rev. Lett. 93 (2004) 012301. 



